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Abstract 

We present results for the masses of pseudoscalar mesons in twisted mass lattice 
, QCD with a degenerate doublet of u and d quarks and a non-degenerate doublet of 

s and c quarks in the framework of next-to-leading order chiral perturbation theory, 
O ■ including lattice effects up to 0{a'^). The masses depend on the two twist angles for 

the light and heavy sectors. For maximal twist in both sectors, 0(a)-improvement 
is explicitly exhibited. The mixing of flavour-neutral mesons is also discussed, and 
results in the literature for the case of degenerate s and c quarks are corrected. 
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1 Introduction 

For non-perturbative studies of Quantum Chromodynamics on a space-time lattice the 
formulation with a chirally rotated mass term [H E] has proven to be a very effective 
framework. This so-called twisted mass lattice QCD implies automatic 0{a) Symanzik- 
improvement if the twist angle is set to a value of tt/2 [21 H]- Numerical simulations 
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of twisted mass lattice QCD have allowed to determine a number of physically relev- 
ant hadronic parameters, notably the Gasser-Leutwyler low-energy constants, to a high 
precision. 

First studies of twisted mass lattice QCD have been made with two mass-degenerate 
quark flavours, representing the up and down quarks, and have allowed to obtain results 
for quantities in the pion sector of QCD [3 El Ej, for a review see In recent work 
these calculations have been extended to 2 + 1 + 1 quark flavours, meaning a degenerate 
doublet of up and down quarks and a non-degenerate doublet of charm and strange quarks 
[SI Uni [m IE] • The chlral twist is implemented with two independent twist angles in the 
up-down sector and in the strange-charm sector. Both of these angles are tuned to 7r/2. 
In the numerical simulations the mass of the charm quark is on the cut-off scale, whereas 
the mass of the strange quark is near its physical value, so that physical results about 
observables in the sector of the the three lightest quarks can be obtained. 

For the analysis of the numerical data, which are obtained at varying quark masses, 
chiral perturbation theory is an invaluable tool, for a review see [13] . Chiral perturbation 
yields analytical formulae for the dependence of physical quantities like meson masses and 
decay constants on the quark masses. It amounts to an expansion around the chirally 
symmetric limit of vanishing quark masses, and is reliable for sufficiently small quark 
masses. 

For the case of twisted mass lattice QCD, chiral perturbation theory has first been 
developed and used for A^f = 2 quark flavours [HJ [T51 [ISl [IZ] . It has been extended to 
A^f = 3, including the strange quark, in [18j. In view of the numerical work on A'^f = 2+1 + 1 
twisted mass lattice QCD, it is the purpose of this article to present chiral perturbation 
theory for this situation. 

Of course, at this point the question of the region of applicability of chiral perturbation 
arises. Chiral perturbation theory is well applicable to quark masses as light as those of 
the up and down quarks. The mass of the strange quark has turned out to at the border 
of the region of validity of chiral perturbation theory, depending on the order of the 
expansion and the use of partial resummations. The charm quark is definitely too heavy 
to be treated within chiral perturbation theory. 

Therefore, chiral perturbation for twisted mass lattice QCD with A^f = 2 + 1 + 1 quark 
flavours is not applicable to quantities in all 4 quark sectors. In this article we restrict 
ourselves to quantities in the pion- and the kaon-sector. Moreover, for physical values 
of the strange quark mass, the formulae cannot be applied to obtain precise quantitat- 
ive results in the kaon sector. Instead, the meaning of these calculations is different. 
Primarily, the idea is to explicitly reveal the structure of A^f = 2 + 1 + 1 twisted mass 
lattice QCD in the sense of its dependence on the two independent twist angles. Chiral 
perturbation displays how the values of the two twist angles affect the observables in the 
different mesonic sectors, and explicitly exhibits the requirements for automatic 0{a) im- 
provement. Secondly, chiral perturbation theory yields the formulae that would allow to 



2 



extract A^f = 3 low-energy constants from data in the region of unphysically small strange 
quark masses. 

In this letter we present the results of chiral perturbation theory for A^f = 2 + 1 + 1 
twisted mass lattice QCD in next-to-leading order. Lattice effects are included up to 
order a?. The special case A^f = 2 + 2 with two degenerate quark doublets has already 
been considered in [12]. In next-to- leading order our results show differences to ^2], see 
below. 



2 Structure of the theory 
2.1 Field definitions 

Chiral perturbation theory is formulated in terms of the pseudo-Goldstone fields for spon- 
taneously broken SU(A^f)j;, ® SU(A^f)i^ chiral symmetry. For A^f = 4, the 15 canonical fields 
0a(2^), CL = 1 ... 15, corresponding to the generators Xa of SU(4) (see App. |A]), are as- 
sembled in the matrix valued field \a4>a{x). Referring to the usual particle multiplets, the 
canonical fields are also denoted {(pa) = (tt, i^, 773, Z), Dg, 7715). The fields corresponding to 
the particle eigenstates are given by the linear combinations 



D± = ^(Dn±i/^i2), 



TT 



D°,5°=^(D9±iDio) 



The field configuration matrix then reads 



V2 



TT 

K- 



TC 



K+ D 

-0 

D- 

Dt 



V2 



[^m + Ai5?7i5] • 



(2) 



It enters chiral perturbation theory in terms of the exponential parameterisation U{x) 



exp 



^ \a(pa{x) 



where Fq is a low-energy constant. The effective Lagrangian for U{x) 
contains chiral symmetry breaking quark mass terms and lattice terms. To leading order 
(LO) it is given by 



£lo = ^ Tr [d,U d,U^) - ^ Tr [xU^ + Ux^) -^Tt[pU^ + Up^) . (3) 



Here the quark mass term is 



X = 25oM 



(4) 
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with the quark mass matrix M, and the lattice artifacts are represented by 

p = 2Woal , (5) 

where a is the lattice spacing and Bq and Wq are additional low-energy parameters. In 
next-to- leading order (NLO) the contributions to the effective Lagrangian contain the 
Gasser-Leutwyler coefficients Li, as well as new coefficients Wi,W-, parameterising the 
lattice artifacts. We refer the reader to pTl [TH] for definitions and details about the NLO 
Lagrangian. 

2.2 Mass terms 

In the lattice regularisation, the untwisted physical quark mass m' is proportional to the 
difference mo — moc of the bare and the critical mass; the prime indicates that the quark 
mass shift already includes an 0{a) correction aWo/Bo, see [T7]. In the (m, (i)-sector the 
twisted mass fi is conveniently introduced via the term ifir^. In chiral perturbation theory 
we use the variables 

X' = 2Bo{m' + i/iTg) = 2Bomg e'"«"^ , = ^m'^ + p? , (6) 

and 

x'o = 25om', X3 = 2So/i, p = 2Woa . (7) 
In leading order (LO) the pion mass is then given by 

ml = 1x1 = \lx'i + Xl = '2Bomg . (8) 

A second twist can be introduced in the (c, s)-doublet in the same way. Considering 
the definition of the fields in Eq. it should be noted that the twist in the (c, s)- 
doublet has the opposite sign, because s and c are exchanged so that they are ordered 
by increasing quark masses in the mass matrix. Following the general argument given 
in [51 S], (9(a)-improvement will be realised at twofold maximum twist, i.e. when both 
twist angles equal 7r/2. Maintaining positivity of the fermion determinant requires to 
implement the twist in the heavy quark sector orthogonal to the mass splitting between 
the heavy quarks [H EO]- Because it is more natural to have a diagonal mass matrix, as 
e.g. in [in], the twist is implemented non- diagonally with ti. 

Referring to the definitions above we define for the light and the heavy sector 

x'o,i = 2Bom[ , X3 = 2BoPi , (9) 

Xo,h = '^Bom'f^ , Xi = 2Boph, xs = Bo{mc - rus) , (10) 

where m'^ is the average heavy quark mass. The symmetry breaking term x' can be 
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written separately for the sectors. With the corresponding Pauh matrices this reads 



, ■ Xo,i +1X3T-3 m {u,d), 

Xo,h + mn + X5T3 in(c, s). 

The LO twist angles are uq^i and uq^h, respectively; for brevity we define cq,; = cos(u;o,05 
co,/i = cos(a;o,/i) and sq,/ = sin(a;o,0) so,h = sin(a;o,h) as in [19]. 

In the following, perturbative quantities denoted by lower case letters (e.g. m^) are 
LO and those with capitals (e.g. M^) are NLO. 



3 Meson masses in leading order 

In the framework defined before, we have calculated meson masses and decay constants 
in NLO chiral perturbation theory, including lattice artifacts up to (9(a^). In this article 
we shall mainly concentrate on meson masses, which are of primary interest. Results on 
the decay constants can be found in f2T\ . 



3.1 Meson masses 

For the meson masses in LO we find the extended Gell-Mann-Okubo mass formulae 



Am 



K 



D 



3m 

6m, 



2 
2 

J?! 5 



+ m^. 



+ m 



D 



m 



(12) 
(13) 
(14) 



where the canonical states r^g and 7715 are not mass eigenstates of the theory. For the 
flavour charged mesons we find in terms of the underlying parameters 

^1 = \/xo' + Xl , 
^l,K = I {^X'o'i + xl + ^X'ol + xl ± Xs) , (15) 

= ^x'o,h + Xi- 

In Monte Carlo calculations, the renormalised physical PCAC quark masses and the bare 
untwisted quark masses can be obtained from the physical and untwisted lattice currents, 
respectively, see [TU]. With their appropriate renormalisation constants on the lattice 
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these full non-perturbative masses are related to each other by [TU] 



mr^^ = ZpV(Z^mPCAC)2 + ^2^ (16) 
±Z-s^^^,. (17) 



m 



PCAC _i_ v-i 



In terms of the PCAC quark masses the perturbative Eq. (I15p reads 



,2 _ r, D ™PCAC 



ml,K = Bomr^^ + B,ml^,^^ , (18) 
3.2 Mixing of flavour neutral mesons 

In general there is mixing between the four flavour-neutral mesons vra, 773 and 7715 from the 
15-plet and the SU(4)-singlet r^i, resulting in the mass eigenstates 7r°, 77, rjc and rj' . Due 
to isospin symmetry in our case, tts is not mixing with the ?7's. In a region with "light" s 
and c quarks, where chiral perturbation theory holds, the rj' can be integrated out and is 
not included in ordinary chiral perturbation theory. 

We define the states 77, rjc in the mass diagonal basis through 



cos 6 sin0\ / r/g 
-sin 6* cos 6*/ 1 7715 



(19) 



The mixing angle ^ is in LO given by 



2 ^2 \ -1 



tan2g = V8 1 + ^ "^; "^f • (20) 
V 2 mt - mi 



The masses of the states in the new basis are then given by 



,2 2 



where tti^^ and tti^^^ are given by the Gell-Mann-Okubo mass formulae ffT^ . ffT^ . and 



(21) 



A.,c = K, - O ^ = - "^') ' (22) 



which is the analogon to the pion mass splitting in completely non-degenerate 3-fiavour 
QCD 
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4 Results in next-to-leading order 



4.1 Vacuum 

In chiral perturbation theory for twisted mass lattice QCD, the vacuum configuration 
is not given by vanishing fields (pa-, corresponding to U{x) = 1; instead there is a non- 
trivial vacuum configuration Uq, which is obtained by minimising the effective action. The 
physical fields represent the deviation of U{x) from the vacuum Uq. It is advantageous to 
parameterise the physical fields through U{x) = y/Uo U-phys{x) \/Uq. In order to calculate 
masses and decay constants it is necessary to determine the vacuum Uq. In our case one 
has 



Uq = exp 



— (a); A3 + cj/iAia) 
^0 



(23) 



In NLO we explicitly find 



m„2W 



nit 



m 



2We] + pco,i m' + pco,/, 4[L6 -W^ + W[ 



(24) 



8 pso,h 



^0 ^Ds 



(25) 

where W = We + W^j^ - 2Le - Lg and W = W;^ + - We - Ws/2 + Lg + Ls/2 

(see m- 

4.2 Pion masses 

Calculating the tree-level and one-loop contributions to NLO, we find for the charged 
pions 

= ml + A(mt[2(L8 + 21^) - (L, + 2L4)] + 2m>^j2L6 - 

- pmlcoMLs + 2Le) - {L, + 2L,) - 2{Ws + 2We) + {W, + 2W4)] 
+ pmlco^h2[L4-Wi] 

- p {m\)CQ^i + mlcQ^h)2[2Le - Wq\ 



+ p'claliLs + 2L6) - (H^8 + 2We) + {W'^ + 2W'^) 
+ p'co,iCo,,,4[L6 -We + W'^]^ + loop,. 



(26) 



The loop contribution is 



loop. 



A/'f(47rFo)2 



(27) 
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with A = 47rF7r conventionally. It depends on 



^^ = ^^{9) = V2sm29-sm'^9 e [0,1] 



(2^ 



The expression for explicitly verifies that (9(a)-improvement of the pion mass requires 
both twist angles to be at maximal twist [21 II]. In particular, a non-maximal twist angle in 
the heavy sector would contaminate the pions with (9(a)-terms proportional to pm'^CQ^h- 
A non- vanishing twist breaks isospin symmetry and leads to a mass splitting between 
charged and neutral pions. For the mass splitting we find 



Ml - Ml = [(Lg + 2L,) - {Ws + 2W,) + (IV^ + 2iy^ 



4.3 Kaon masses 

The kaons are degenerate at NLO. For their masses we obtain 
4 



(29) 



Mi 



m 



K + ^(2mi[2(L8 + 4L6) - (L5 + AL^)] + Am\{ml - m\)[2L^ - L,] 

- pml{co,i + Co,hmLs + 4L6) - {L, + 4L4) - 2{Ws + AWf,) + {W^ + AW^)] 

- p{m%~ mj^){co,i + Co,h)2[2LQ - Wq] 
+ p\co,i + co^hfliLs + ALe) - {Ws + AWe) + {W^ + AW^,)] 

- P'(4/ + ^o,h)[^8 - + W^g]} + looPi. (30) 



with 



and 



loop 



K 



K 



NiiA-KFof 



4-6 



mz 



~ 62 2 1 

m„ In 

3 



mz 



A2 



^2 = 62(^) =5sin2^- v^sin^cos^ G 



5-3^3 



(31) 



(32) 



Again, 0{a) improvement is obtained when both light and heavy twists are maximal. In 
contrast to the case of pions, however, the kaon masses contain 0{a?) terms at maximal 
twist. 



4.4 Heavy-light case 

The heavy-light case is the situation where the both sectors are mass degenerate, i. e. 
mi = mu = md and mh = = mc- In this special case the mixing angle is analytically 
given by 

tan 20 = 78 (33) 

This mixing angle corresponds to including the physical states rj and rjc in the field matrix 
by replacing the generators Ag, A15 by Ag, A^g (see App. |A|) . 
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In LO the mass spectrum given in Eqs. fllSp . 0181) . (^T^ simplifies to 



K 



m 



D 



m 



Ds 



PCAC 



25omf 



(34) 



The loop contributions simplify too, because of 



(35) 



In NLO there are mixing effects between the K- and D-mesons, when = rric- In the 
generic case 7^ rric these mixings contribute only at NNLO, related to the non-analytic 
behaviour in — rric. In the case = rric, however, the mixing already shows up in 
NLO. In particular it leads to a fiavour breaking in both sectors according to 



M|± - M|.o = Ml± - Mlo 



(36) 



and the above formula fl5Ul) for M^- gives the average of the kaon masses. This breaking 
is a pure lattice artifact. 

At tree level the resulting formulae for the kaon masses coincide exactly with those 
given in [12]. In NLO, however, the loop contributions to the kaon masses 



A^f(47rFo)2 



mjlnl^ 



(37) 



which we obtain, are different from those given in [19j : 



m 



K 



iVf(47rFo)2 



-m^. In 



m. 



A2 



1 
3 



m~_ In 



m. 



V15 



A2 



(3^ 



The origin of this discrepancy is the neglect of the mixing of the 77-mesons in [12]. The 
results are only valid if the mixing disappears, which only applies if all quark masses 
are degenerate. Also, as a consequence, the loop-divergences in [TH] do not cancel the 
tree-level ones if mi 7^ m^. 

A related discrepancy can be found for the kaon decay constant, where our result for 
the loop contribution is 



2(47rFo)2 



3 21 /< 
— mi In — ^ 
16 ^ V A2 



+ -mj^ In 



m 



K 



A2 



+ —m^^ In 



m 



16 



A2 



(39) 
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in contrast to P^ : 



2(47rFo)2 



+ In 



m 



K 



A2 



m: 



A2 



m 



Ds 

A2 



(40) 



5 Conclusion 

We have obtained the masses of pseudoscalar mesons in the framework of chiral perturb- 
ation theory for twisted mass lattice QCD with a degenerate doublet of u and d quarks 
and a non-degenerate doublet of s and c quarks in next-to-leading order, including lat- 
tice effects up to C(a^). The results display the dependence of the masses on the two 
twist angles for the light and heavy sectors. For maximal twist in both sectors, Symanzik 
improvement to 0{a), as proven in [2111], is explicitly exhibited. The mixing of flavour- 
neutral mesons is also discussed. For the case of degenerate s and c quarks, proper account 
of mixing corrects results in the literature. 



A Generators of SU(iV) 

The generators of SU(A^) in the fundamental representation are the hermitian N x N 
matrices with zero trace. A useful representation is as follows [23] ■ Let Cij be the matrix 
with the only non- vanishing entry being 1 in row/column ij: 

{Gij)ki = SikSji . (41) 

The off-diagonal generators of SU(A^) can be written in terms oi N x N generalisations 
of the Pauli matrices ai and (T2, defined by 

cr2;ij = -^eij + icji , {i < j). (42) 

The diagonal traceless matrices, generalising a^, are defined by 



^3;^ = y ^(-^ diag(l__J., 1 - z, 0^ , 0) , (z > 1). (43) 

i~l N- 

The extended Gell-Mann matrices are then given by 



(44) 





= 5 


for 


a 


= (J- 


l)2 + 2i-2. 


< j), 




= Aa , 


for 


a 


= (J- 


lf + 2t-l, 




C^3;i 


= , 


for 


a 


= z'- 


1, {t>l), 
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where a = l,...,iV^ — 1, and the generators of SU(A^) are equal to Aa/2. They are 
orthonormal in the sense of 

Tr(AA) = 2(5,,. (45) 

For SU(4) the diagonal Gell-Mann matrices, corresponding to the flavour neutral 
mesons, are explicitly given by 

As = diag(l, -1, 0,0), 73 As = diag(l, 1,-2,0), Vq A15 = diag(l, 1, 1, -3) . (46) 

In the heavy-light case the latter two can be replaced for simplicity with the linear 
combinations 

V2Xl = diag(l, 1, -1, -1) , XI, = diag(0, 0, 1, -1) . (47) 
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